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Abstract 

We undertake first steps in making a class of discrete models of quantum gravity, spin 
foams, accessible to a large scale analysis by numerical and computational methods. In 
particular, we apply Migdal-Kadanoff and Tensor Network Renormalization schemes to spin 
net and spin foam models based on finite Abelian groups and introduce 'cutoff models' to 
probe the fate of gauge symmetries under various such approximated renormalization group 
flows. For the Tensor Network Renormalization analysis, a new Gaufi constraint preserving 
algorithm is introduced to improve numerical stability and aid physical interpretation. We 
also describe the fixed point structure and establish an equivalence of certain models. 
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1 Introduction 

Spin foam models aim at providing a description of the microscopic structure of spacetime and 
thus a theory of quantum gravity [1-11]. These models can be understood as a non-perturbative 
definition of the path integral for quantum gravity. To make these path integrals well defined 
one has to introduce a regularization based on a choice of discretization, i.e. a lattice or, more 
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generally, a triangulation or two-complex. Indeed, spin foams can be understood as generalized 
lattice gauge theories. 

This discrete (a priory auxiliary) structure should not be confused with another feature of 
spin foams, which is often termed 'Planck scale discreteness' [12-16], namely that the spectra of 
(kinematical) geometrical quantum observables, like areas and volumes, are discrete. There are 
thus two different kinds of UV cutoffs, whose interplay has not been fully understood yet. This 
has to be kept in mind when discussing a possible breaking of (global) Lorentz or (local) diffeo- 
morphism symmetry. A (naive) lattice regularization will generically break these symmetries, 
see for instance [17-24] for a discussion of these issues in gravity. 

We may, however, consider a continuum limit with respect to this auxiliary discretization 
scale, for example by a coarse graining or blocking procedure, see [25-27] for recent examples 
involving gravity or related to it. A crucial question then is whether Lorentz or diffeomorphism 
symmetry will be restored in this limit, despite the possibility of still having the second kind 
of UV cutoff, provided by the discreteness of the spectra, in the theory. That this cutoff does 
not necessarily lead to a violation of Lorentz symmetry has, for instance, been argued in [28] 
on kinematical grounds. A full dynamical scenario for 4D gravity where a restoration has been 
shown to occur is, however, missing, nonetheless see [29, 30] for progress in this direction.^ 

These questions motivate us to consider a continuum limit which involves many building 
blocks or large lattices (with many vertices), as this is the limit where one can hope to obtain 
a diffeomorphism invariant theory. An alternative is the semi-classical limit [36-38], in which 
rather the Planck constant, leading to 'Planck scale discreteness', is taken to zero. To distinguish 
these two kind of limits we will sometimes refer to the first one as statistical limit. 

Experience with other quantum gravity models, such as (causal) dynamical triangulations, 
has shown [39-45] that even before the question of restoration of symmetries can be addressed, 
it is not at all obvious whether such a statistical limit leads to any viable model of spacetime, i.e. 
whether such a limit will result in smooth four-dimensional spacetime manifolds. Indeed, in this 
kind of limit statistical considerations become important and it can easily happen that state sums 
become entropically dominated by configurations not resembling any four-dimensional manifold 
at all. As we will see, a related issue arises for spin foam models (or other models based on first 
order /tetrad formulations) where geometrically degenerate configurations might be dominant. 
Such configurations also turn up in a semi-classical or classical phase space analysis, even if this 
involves only a single simplex [46-48]. 

Hence it is crucial to investigate which kind of large scale physics or, in other words, phases, 
are encoded in the candidate quantum gravity models. Phases are often characterized by sym- 
metries, that is such a study might also answer which kind of symmetries might be restored 
in a large scale / statistical limit. Indeed, making progress in this direction is one of the most 
pressing issues for the spin foam approach. However, it is also a long standing open issue [49-51] 
charged with a number of conceptual and technical challenges. 

One important challenge is the complex structure of the models which lead to very com- 
plicated amplitudes as compared, for instance, to QCD. Here, our strategy [11] is to develop a 
wide range of simplified models which capture essential features of spin foams while being much 
easier to handle. These simplifications are obtained on the one hand by replacing Lie groups, 
on which the gravitational spin foams are based, with finite groups. On the other hand, we can 
also consider 'dimensionally reduced' models (spin net models). In fact, 2D spin net models of 
the simplest class share many statistical properties with their counterpart 4D spin foam models. 

Similar simplified models have been successfully studied, e.g. [52-54], to get insights into the 

^For 3D gravity, which is a topological theory, that is without propagating degrees of freedom, discretization 
does not necessarily lead to a breaking of diffeomorphism symmetry [20, 31, 32]. This holds also for 3D gravity 
with cosmological constant [33]. One can, however, consider discretization or quantization methods which a priori 
break diffeomorphism symmetry and look for methods to restore these symmetries, see [25, 34, 35]. 
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large scale behavior of lattice gauge theories. Here we hope for a similar improved understanding 
of the possible phases that can occur in quantum gravity models. In particular, we will see in the 
course of the paper, how conjectures or even conclusions for Lie groups can be made based on 
findings for finite groups. These simplified models can also be interesting in their own right [55- 
57], in particular if an example is found in which some analogue to diffeomorphism symmetry 
is restored. Indeed topological phases and string net condensates [58] which are studied in 
condensed matter, also regarding the question of symmetry restoration, are tightly related to 
3D spin foams with finite groups [11]. 

The development of coarse graining and renormalization techniques seems to be the most 
promising avenue to study the large scale behavior and simplified models allow us to adapt and 
further refine methods from lattice gauge theory and condensed matter systems. In this work we 
will therefore apply the Migdal-Kadanoff scheme [59, 60] and the tensor network renormalization 
(TNR) method [61, 62]. These schemes involve a regular lattice and, due to this regular structure, 
they are amenable to efficient numerical simulations. 

With this approach we are able to explicitely answer the question of BF symmetry restora- 
tion for a range of models and to gain insights into how these results are related to the case of 
infinite groups. These results should be understood as a first step towards harnessing the power 
of numerical methods from statistical physics to deepen the understanding of the large-scale 
physics of spin foam models. 

In this work, we will first introduce spin foam and spin net models and write them in ways 
suitable for coarse graining (section 2). We will also define a particularly important class of 
models, termed 'Abelian cutoff models', discuss the role of BF / translation symmetry and 
detail the relationsship between spin foams and nets. Subsequently, in section 3, we discuss 
the conceptual challenges of coarse graining in this context and argue for the approach pursued 
here, in particular for the use of a regular lattice. 

We then apply the Migdal-Kadanoff and tensor network approximation schemes to coarse 
grain our models (sections 4 and 5, resp.). In each case, we first introduce the method and 
highlight some of its analytical properties before presenting numerical results that focus on 
the question whether renormalization of Abelian cutoff models will restore BF symmetry. In 
particular, we feature a Gauf5 constraint preserving TNR algorithm tailored to the geometric 
interpretation of our models and establish an equivalence amongst certain cutoff models under 
the TNR renormalization scheme. We conclude by comparing both approximation schemes and 
pointing out possible future directions of research. 

2 Spin foams and spin net models 

Spin foams are a particular class of lattice gauge models (see e.g. [63] for a recent review and 
[11] for a review emphasizing the relation to lattice gauge and statistical physics models). Such 
models are specified by variables, taking values in some group G, associated to the edges of 
a lattice (or more generally an oriented 2-complex) and weights associated to the plaquettes. 
They can thus also be termed plaquette models. 

A related class of models, which will be introduced below, are so called edge or spin net 
models [11]. Here group variables are associated to the vertices of a lattice (or more generically 
an oriented graph or 1-complex) and weights to the edges. This class includes the well-known 
Ising models, based on the group Z2. Indeed it will turn out that the structures involved in a 
spin net model are very similar to those involved in spin foam models - just that where, for 
instance, weights are associated to 2D plaquettes for spin foams, weights are associated to ID 
edges in spin nets, similarly for the group variables and so on. In this sense spin nets are a 
simpler or dimensionally reduced form of spin foams. There is however one essential difference, 
which is that spin foams enjoy a local gauge symmetry whereas spin nets only feature a global 
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symmetry, in both cases given by the group G the models are based on. In section 2.3 we will 
also comment on another relationship between spin foams and spin nets: spin nets can be seen 
as measuring non gauge-invariant observables in spin foam models. 

Spin foams and spin nets are defined by partition functions, and we will first consider a 
representation of these partition functions as sums over group variables. Via a group Fourier 
transform we can rewrite these partition functions as sums over variables labeling the irreducible 
representations of the group G. This is where the name 'spin foam' stems from, as 'spin' refers 
to the representation labels for the group SU{2). This alternative representation is well known 
as a duality transformation for both edge and plaquette models, and is usually employed for 
the high temperature or strong coupling expansion [64]. The models in this representation 
are not only specified by the dual weights but also by an intertwining projector acting on a 
certain representation space. Non-trivial spin foam or spin net models can be constructed by 
choosing this projector to be different from its standard form (which is the Haar intertwiner 
introduced below) in plaquette and edge models respectively. In the case of Abelian groups 
we will explicitly construct non-trivial models, the so-called Abelian cutoff models, which only 
consider representation labels of the group up to a certain cutoff. For these models, we will see 
that the choice of a non-trivial projector is equivalent to retaining the Haar intertwiner while 
restricting the dual weights in a particular way. These are precisely the models that we will 
numerically analyze in sections 4 and 5. 

For the non-trivial models, one can then apply the inverse group Fourier transform and again 
obtain a partition function in terms of group variables. As will be explained below, this will 
however require the introduction of several group variables per vertex (for edge models) or edge 
(for plaquette models) [65-67] . This representation is termed holonomy representation, both for 
spin foams and spin nets. 

In the next two subsections we will give a short introduction into the main concepts and 
different representations of spin foams and spin nets. Furthermore we will detail the different 
possibilities of rewriting these models into tensor network form, as this will be the basis of one 
of our coarse graining methods, to be discussed later on. 

2.1 Spin net models 

To construct spin net models we start with state sum models formulated over a finite group 
G and on a graph (one-dimensional complex) with oriented edges. More precisely we consider 
partition functions of the type 

' ' {9v} e 

where (ju is the number of vertices in the graph, s(e) denotes the source vertex (starting point) 
of the edge e, t{e) denotes its target vertex (final point), and the curly brackets under the sum 
symbol denote that there is a sum per vertex: X^j^^} = HoX^g^- Here group elements are 
associated to vertices and weights, which determine the couplings, to edges. Therefore these 
models are also known as edge models and include the standard Ising model for which G is 
equal to Z2. The weights We{gs{e)9^(l)) arbitrary functions^ over the group G. However, 

if We are class functions, i.e. invariant under conjugation {we{g) = We{hgh^^) yg,h £ G), the 
model will feature a global symmetry given by the group G: the partition function remains 
invariant when applying the same conjugation to group variables at each vertex. 

The form (2.1) defines the simplest form of spin net models in the representation based 
on group variables. This representation will be called holonomy representation in analogy to 

^To have a statistical interpretation of We as probability weights these should be positive. This will however 
not necessary be the case for spin foam models. 
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Pe{9v)a,c, e P*e{9v')b,c, Pe{9v) P*e{9v') 

Figure 1: On the left: the three objects associated to every edge. On the right: their schematic 
representation. Every straight hne joining two objects means a contraction of indices. 



spin foam models (where group variables are associated to edges and represented holonomies 
of a connection). Holonomy representations of more general models will require several group 
variables associated to each vertex, as we will see later. 

Via the group Fourier transform, we can change from the above representation in terms of 
group elements, to the spin net representation, which is in terms of the irreducible representa- 
tions of the group^ . Every function on the group can be decomposed in terms of matrix elements 
of the irreducible representations p, 

dimp ^. 
^(9) =J2Y1 {Wp)abP{9)ab, {Wp)ab = "T^ ^ w{g) p*{9)ab, (2.2) 

being p* the dual of p. For class functions this decomposition reduces to the character decom- 
position 

w{9) = ^WpXp{9), Wp = —^w{g)xp''{g), (2.3) 
p g 

where Xp{9) = Yl^=i P{9)aa denotes the character. We note that our convention for the delta 
function over the group, Sg, is 



^ E ^G{g)f{g) = fi^d), 5g{9) = E dim/9 Xp(5). (2.4) 
g p 

Using the property 

P{9'^)ab = P*{9)ba. (2.5) 

we obtain ^ 

Z = X] n (*Pe)aefee Pe{9s{e))a,c, P*e{9t{e))b,c,, (2.6) 

{gv}{Pe} e 

where we sum over repeated indices.^ Note that, associated to every edge, there is a coefficient 
{wp^)aeb^, and two group representations, Pe{9s{e))a,^c,^ living on the source vertex and pl{gt(e))beCe 
living in the target vertex. The indices of these three objects are contracted, as described 
schematically in figure 1. 

Now we can carry out the sums over group variables. The result is 

^ = E (n«)'^^^e) n^ae...,V...;Ce...,c.,...({/'e}e^.) (2.7) 

{pe} e 

where we have defined the vertex weight 

Pl...,b^,...;c,...,c^,...i{pe}eDv) ■■= J^^Yl P^(9v)a,c, P*e'^9v\,c^, ■ (2.8) 

gv e:v=s{e) e':v=t{e') 



^We refer the reader e.g. to [68] for the main concepts of group representation theory that we employ. 
*For class functions We ^ 
(pe)aeCc to the characters Xp. 



*For class functions We we have (wp^)a^b^ = (w)p^5a^b^, which will contract the representation matrices 
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Pe2 {9v) 



K P*eA9v) 



P*eA9v) 



^64 *^e4 



P''(Pei,-,Pe4) 



^64 ^^64 



Figure 2: Four-valent vertex with two outgoing edges, ei and 62, and with two incoming edges, 
63 and 64. On the left: representations meeting in the vertex. On the right: schematic repre- 
sentation of the resulting vertex weight. 



The first and third groups of indices involve all the edges for which v is the source vertex. The 

second and fourth groups of indices involve all the edges for which v is the target vertex. Using 
the schematic representation employed in figure 1, we can represent the vertex weight as in 
figure 2. 

Note that P'" can be seen as an intertwiner map, called the Haar intertwiner, acting on a 

certain representation space for the group G. This representation space, Ti^, associated to the 
vertex v, is given by the tensor product of the representations pe associated to the outgoing 
edges and the representations p*, associated to the incoming edges 




The intertwining property of this map is guaranteed by the averaging over the group in (2.8). 
Indeed, the Haar intertwiner defines an orthogonal projector onto the subspace W^"" of Hv, 
invariant under the group action defined on this representation space. 

The same intertwining map will appear in spin foam models. There, the choice of the Haar 
intertwiner as a projector and face weights to be trivial defines topological models, known as 
BF-theories. A gauge symmetry, known as translation symmetry, arises in this case, forcing 
the model not to have local physical degrees of freedom. The analog situation happens in spin 
net models. The choice of edge weights We = id and of the Haar intertwiner as a projector 
corresponds to the model at zero temperature, with no local degrees of freedom. For the weights 
We this amounts to choosing We ^ 6g- In this case the projectors are just contracted along 
the edges of the graph. 

Non-trivial models are more interesting. They are constructed by restricting the projector 
further, i.e. by selecting a subspace of the invariant subspace of Hy and by replacing the Haar 
intertwiner (2.8) by a projector onto this subspace. We will proceed in that way here, and in 
general assume that P^ is a projector onto some subspace of 7{*"^. This allows us to obtain 
interesting models even if we choose the edge weights We to be trivial. Indeed, in these models, 
the original gauge symmetry is broken and local physical degrees of freedom arise. We will see 
this behavior for the Abelian cutoff models described below. 

In general, in order to re-express the partition function Z of such non-trivial models in the 
holonomy representation, namely as a sum over group variables, we will need to associate more 
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than one group elements to each of the vertices in the lattice. Let us denote the number of edges 
attached to the vertex v by n^. Now, we assign one group element g{v, e) to each of the edges 
attached to v and define 

^''({5(.,e)}e3.) = E^ae...,6.,..,c....,c,,...({/^e}e3.) 

X Yl ^^Mpe) P*e{g{v,e))a,c, '^^Mp'e) Pe'{g(v,e')\,c^, ■ (2-10) 

e\v=s{e) e'\v=t{e') 

In terms of this vertex amplitude, the partition function reads 

^=(11^) E (n^^(^We).e)5(-;fe),e)))(n^^(i5Ke)}«^))- (^-H) 

^ {!?(.,e)} e V 

In case that is given by (2.8), i.e. by the Haar intertwiner, we obtain the form of the partition 
function given in (2.1), as P'" then enforces equality between the group elements associated to 
one and the same vertex. 

A different simplification occurs when We{h) ~ dcih), a choice that we already mentioned 
above. Then the two group elements g(v,e) g(^v',e) associated to any edge e have to be equal. 
Hence, the partition function reduces to a sum over group elements ge = 9{v,e) = g{v',e) associated 
to edges. This type of models is known as vertex models - the energy of a configuration is now 
determined by the vertex weights P^{{ge}eDv)- 

In this work we will apply coarse graining to models with Abelian groups Zg. In this case, 
all irreducible unitary representations, which we will label by A; G Zg, are one-dimensional and 
defined by their characters Xk{g) = exp(^/c • g) for g G Zg. The transformation between 
functions on the group Zg and on the dual, equal to the space of characters, which is also given 
by Zq, is given by the discrete Fourier transform 

q-l q-1 

w{g) = ^WkXk{9) Wk = q~^^w{g)xk{9)- (2.12) 

A:=0 9=0 

Characters for Abelian groups are multiplicative, i.e. Xk{9i '92) = Xkidi) " Xki92), and also 
Xk{g^^) = Xk^ig) — Xkig)- Moreover, the delta over the group is now the g-periodic delta. It 
verifies q~'EgS^'H9)f{9) = fiO). 

The spin net representation of the partition function simplifies in the case of Abelian groups 

to 

{ke} e ^ <^'^^ 63*^ 

where is equal to +1 (—1) if v is the source (target) of e. The projector P'" implements the 
GauB constraints at the vertex v. It indeed projects to the irreducible subrepresentation in the 
tensor product of all representations associated to the outgoing edges and the tensor product 
of dual representations of incoming edges. The tensor product is one-dimensional and equal to 
the trivial representation if the oriented sum of the representation labels fee is equal to zero. 

Note that the spin net representation is the starting point for the high temperature expansion 
[69]. The infinite or high temperature fixed point is represented by We{k) = 6^'^\k). On the 
other hand the zero or low temperature fixed point is given by We{k) = 1. This corresponds to 
weights We{h) = dcih) in the original group representation. The model is 'frozen', i.e. the group 
elements at the different vertices have to agree (assuming that the graph has only one connected 
component). 
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As commented before, for this choice of weights there is a gauge symmetry. It is associated 
to the faces, i.e. the two-dimensional cehs (here we are assuming that the graph is actually given 
by an orientable 2~complex, i.e. the 2-dimensional cells are well defined). Associating to every 
face / an element kf G Zg we define a gauge transformation acting as 

ke^k'^ = ke + ^elkf (2.14) 

where e( is +1 (—1) if the orientations of e and / agree (disagree). Under such a transformation 
the contribution of a configuration {ke}e to the partition function Z does not change. Choosing 
either the edge weights We or the vertex projector P'" to be non-trivial will in general break this 
translation symmetry either completely or down to a smaller symmetry. Choosing a non-trivial 
projector will in general result in vertex models, as one can basically reduce the set of vertices 
allowed by the GauB constraints even further. 

We are now in position to introduce the Abelian cutoff models. As said before, for the low 
temperature fixed point (analog to BF in spin foam models) the weights are 

We{9) = S'^'H9) = Y.^>^(9) ^ we{k) = iyk. (2.15) 

{k} 

The Abelian cutoff models are derived from this by 'cutting off' the sum at some value K such 
that some of the dual weights We vanish. Explicitly, 

f 1, for|A;|<K 
Mk) = l ' ^ ', ^ 2.16 
[ 0, for \k\ > K 

where we will consider only even q, hence K < ^. Here the range for k is given by ^ < /c < |. 
Also, note that the symmetry condition w{k) = w{—k)\/k is fulfilled. This requirement is 
desired in the quantum gravity setting because it ensures that the model does not depend on 
edge orientations and certain types of face and edge subdivisions [70, 71]. 

Abelian cutoff models could be equivalently described by a restriction of the projector P^ 
(and keeping the weights We{k) = 1). They provide a simple example of breaking the translation 
gauge symmetry from the frozen model in order to introduce physical degrees of freedom. Also 
this choice of model corresponds to a regularization one often choses for BF lattice theories 
with Lie groups [72] . In this case the orbits of the translation gauge symmetry are non-compact 
and the evaluation of the partition function gives generically infinity. Different methods of 
regularization have been developed, one would be equivalent to introducing a cutoff K (for a 
theory with group G = U{1)). One can then ask whether these regularized models would fiow 
back to the full BF model, or more generally the low temperature fixed point, under coarse 
graining. We will consider this question in sections 4 and 5. 

So far we have presented the partition functions for spin net models as sums over group 
variables (the holonomy representation) or representation labels (the spin net representation). 
An alternative is to write the partition function as a contraction over tensors attached to the 
vertices of the underlying lattice (or some graph associated to the lattice), that is, in the tensor 
network representation. The tensor network representation is commonly employed in statistical 
and quantum systems [61, 73-76], since it is especially suitable for developing techniques of 
renormalization. We will make use of this representation of the partition function in section 5, 
where we will apply the renormalization approach to spin net models with Abelian groups. 

In this representation, the contraction of the indices is prescribed by the edges. The partition 
function can hence be expressed tensor trace 

^ = En^{Me.. = tTr0.r^ (2.17) 

{fee} « 
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Figure 3: Employing the schematic representation of figures 1 and 2 it is straightforward to 
reaUze that the partition function of spin net models can be written in the form of a tensor- 
trace over a network of tensors, being the tensor-trace the sum over representations. 




Figure 4: Schematic definition of the vertex tensor T^^^ | ^ in a square lattice. 



This way of representing the partition function is related to the so-called graphical calculus 
[77-80], which is often employed in the spin foam literature. 

In the particular example of Abelian models in the spin net representation, we first absorb 

~ 1/2 

the edge weights w^^ into the vertex weights -P^({/ce}eDii) by distributing Wj^ factors to each 
of the adjacent vertices, namely in every vertex we define the tensor 

%.U^. = ( U(^kJk'^')P''i{ke}e^.) (2.18) 

The partition function is then given by the tensor trace (2.17) (see figure 3 and figure 4). 

More generally, we can find tensor network representations also for the non-Abelian models 
in the different representations. The holonomy representations is a convenient starting point 
for a low temperature expansion, the spin net representation to the high temperature region. 
For the spin net representation we can proceed as for the Abelian models. For the holonomy 
representation assume that we have edge weights such that We{h) = Scih). We then just 
need to understand the group elements ge as indices attached to the tensors T^^^^ = and 
= to see that (2.11) can be rewritten as a tensor trace. For the more general case of 
non-trivial edge weights we can introduce another set of rank two tensors w to the midpoints 
of the edges. 
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(a) definition of curvature 



(b) edge and faces attached to it 



Figure 5: (a) A face / bound by the edges ei, • • • , 65. The curvature is hf = geige297^ 9eA9e^ ■ 
(b) An edge e with four faces attached, /i and /4 with positive relative orientation with respect 
to e, and /2 and /a with negative relative orientation. The corresponding edge- weight is denoted 

by -f 1 f4 f2 f3 fl fi f 2 f3 ■ 

n-" /T* n-' /!■' •/t' n-' /r* /r' 



2.2 Spin foam models 

Spin foams are in many respects similar to spin nets. The main difference is that spin foams arc 
gauge theories, formulated with a gauge group G, which here will be a finite group. Furthermore 
spin foams require an oriented two-complex for their definition. This implies, that we have a 
well defined notion of oriented edges as well as oriented faces (which are the 2D cells of the 
complex, or the plaqucttcs for a regular lattice). 

In order to define a spin foam model, we assign a group element to every edge e of the 
two-complex and a weight Wf : G ^Cto every face /. The state sum model is then defined by 
the partition function 

^=^En^/(M> (2-19) 

where t|e denotes the total number of edges in the two-complex. 

The function Wf \sa, class function. Furthermore, in (2.19) Wf depends on the group elements 
only through the holonomy hf around the closed loop of edges forming the face /, that we will 
call curvature. Let us make its definition explicit. For that, we recall that the relative orientation 
between a face / and any edge e in the boundary of / is denoted by ei, and it is equal to +1 
(— 1) when e and / have the same (opposite) orientation. Given a face / bounded by the ordered 
sequence of edges ei, 62, • • • , e„, the associated curvature is given by 

hf = gi'gi---gi\ (2.20) 

as depicted in figure 5a. These properties of the weight Wf guarantee that the partition function 

is invariant under gauge transformations ge — )■ gs{e)ge9^(ly A-S before, s(e) denotes the source 
vertex of the edge e while t{e) denotes its target vertex. 

The partition function (2.19) describes standard lattice gauge theories. The weights can for 
instance be chosen to emulate the Wilson action 

wf{h) = eM-Sw{h)) , Sw = ^{Xp{h) + Xp{h-')) , (2.21) 

where a is a coupling constant. For the choice Wf{hf) = daihf), with 6g the delta function 
over the group defined as before, the partition function only sums over (locally) flat holonomies. 
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This is a discretization of BF theory, which is a topological field theory (without propagating 
degrees of freedom). It coincides with the zero temperature fixed point or zero coupling fixed 
point of lattice gauge theories of Yang Mills type. 

In order to obtain a representation of the partition function (2.19) as a sum over represen- 
tation labels (spin foam representation), we again apply the group Fourier transform. Here we 
only need to decompose class functions into characters 

^(5) = '^pxpig), = 1^ H ^(9) Xp'ig) ■ (2.22) 

p ' ' 9 

el el e-f 

To decompose Wf{hf) = Wf{gel^ge2^ • ■ ■ 9eT) use the properties 

Xpigh) = Yp{9)abp{h)ba, P{g~^)ab = P*{g)ba- (2-23) 

ah 

We introduce this decomposition in equation (2.19) and individually carry out the sums over 
the group elements (note that given an edge e, the groups elements ge and g^^ appear in YljWf 
as many times as number of faces share the edge e) . We obtain the following expression for the 
partition function 

dimpy 

^=E(n*p/) E J- v+ (w}/3e), (2.24) 

"s r--!"s i---i"t i---)"t !••• 

{Pf} f {4=l},cf " 

where, associated to every edge, we have defined the projector 

J- ({M/3e):=r^E n n ^/-(5e)„/-a/-- (2-25) 

as ,—,as ,...;a( ,...,a( |(jr| ^ — ' -^-^ «s «t •' a* "s 

In the above tensor, the first and third groups of indices (distinguished by the superindex /"'") 

involve all the faces that have e in their boundary with the same orientation as the face. These 

groups of indices have therefore as many indices as number of faces with el = 1. In turn, the 

second and fourth groups of indices (distinguished by the superindex /~) have as many indices 

as number of faces / with ei = —1. We show an example in figure 5b. To every pair face-edge 

f f 

two indices are associated, ai and ai , that we attach to the vertices of / that bound the edge 
e, namely s(e) and t{e). In equation (2.24), for every face / of the two-complex, there is a sum 
for every vertex v belonging to that face. Note that every index al appears twice, since there 
arc two edges meeting at the vertex v and bounding the face /. Then, the product over edge 
projectors contracts all the indices. 

Associated to every edge there is a Hilbert space Tie- Let us denote by /i, • • • /m the faces 
attached to e. Then, 

He := Vp^^ ®Vp^^®...® Vp^^ . (2.26) 

Here, to keep the notation simple, we are assuming that £6 = 1 for all the faces^. The corre- 
sponding tensor defines an orthogonal projector onto the gauge-invariant subspace of the 
edge-Hilbert space T-Le 

(P^) (2-27) 

^If there would be some face with opposite orientation to that of e, we would replace the corresponding vector 
space Vp. for its dual V* . 
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for -0 G T-Le- Here, as for the spin net models, the projector (2.25) is the Haar-intertwiner. As 
before, more general spin foam models are constructed by restricting the Haar intertwiner to 
proper subspaces of the gauge invariant subspace of He- 

For spin foam models one usually reorganizes the partition function (2.24) such that am- 
plitudes can be associated to vertices. To this end, one decomposed the projectors P'^ by 
introducing an orthonormal basis t'^*, k = 1, . . . ,m, for the invariant subspace of He- By ad- 
justing the basis we can decompose any gauge invariant projector as (here m' < m and m' = m 
for the Haar intertwiner) 

m' m' 

-^'^ = X] I I) -^aiaa-a^,; f)i62-6„/ = ^-a^aa-a^/^'biba-b™^/' (2.28) 

k=l k=l 

SO that the indices attached to t{e) are now carried by the intertwiners and the indices 
attached to s{e) are carried by the intertwiners According to the description above, we can 
contract within every vertex the corresponding obtaining vertex amlitudes Av{pf^L^), that 
depend on the representations pf and intertwiners associated to the faces and edges that meet 
at V. With this, the partition function can be written in terms of vertex amplitudes via 

^ = E n^p/ WMpf.n- (2.29) 

This is the usual description employed in spin foam models for quantum gravity. 

As for the spin net models, we can also define a holonomy representation [65-67] for spin 
foams using the inverse group Fourier transform. In case we are dealing with a non-trivial edge 
projector (2.28), i.e. not the Haar intertwiner, the resulting holonomy representation will be 
of the form of a generalized lattice gauge theory. That is, instead of only one group variable 
associated to every edge we will have as many group variables attached to a given edge as there 
are faces attached to this edge. Furthermore, there will not only be face weights Wf but also 
edge weights P'^ which result from the transformed edge projectors P"^. 

Let us now consider the case of Abelian groups Zg. As for spin net models, the spin foam 
representation simplifies since the irreducible representations are just one-dimensional and the 
group Fourier transform is just given by the usual discrete Fourier transform. Namely, we obtain 

{kf} f e /De 

{kf} f V e^v fDe 

where the g— periodic delta function 6^'^\-) enforces the Gaufi constraints, now based on the edges 
(instead on the vertices as for spin nets). In the last step we just splitted the delta functions over 
every edge into two delta functions over every vertex, in order to define the vertex amplitudes 
(which here are A, = Ue^J^''\Ef:,edkf)). 

In case the face weights are given by = 1, which is the BF theory case, we obtain an 
additional symmetry for the partition function (2.30), which as commented before is known as 
translation symmetry. For spin foams the gauge parameters kc G are based on the 3-cells 
of the lattice (here we assume that the 3-cells are well defined, for a regular hypercubic lattice 
these would be the 3D cubes): 

kf ^ k'f = kf + e)kc (2.31) 
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Figure 6: Schematic definition of the tensor network for spin foams on a cubic lattice. 



where = +1 (= —1) if the orientations of the 3-cell c agrees (disagrees) with the one of the 
2-cell /. As for the spin net models this translation symmetry does not change the validity of 
the Gaui3 constraints appearing in the partition function (2.30). For the gravitational spin foam 
models (in 3D) this type of gauge symmetry gives the diffeomorphism symmetry underlying 
general relativiy [20, 31]. It hence has a special status and is indeed deeply intertwined with 
triangulation or more generally discretization independence of the models [20-22, 25, 27]. 

For coarse graining we will consider the Abelian cutoff models with face weights coinciding 
with the edge weights in spin net models (2.16). For the cutoff models the translation symmetry 
will be broken as the weights are now no longer constant in kf. The motivating question will be 
to see whether these models flow back to the BF phase, in which the translation symmetry is 
restored. 

Finally let us note that also spin foams can be represented in various ways as tensor networks. 
One possibility would be to start with the representation (2.29) involving vertex amplitudes. To 
obtain an algebraically similar form to the tensor network representation of spin net models one 
would however keep the edge projectors P*^ intact. To this end we associate the tensor T*^ = 
to the midpoints of the edges of the lattice. The edges of the lattice carry a number of indices 
which are all contracted with each other in the lattice vertices according to the description below 
(2.25). Furthermore we have to take care of the face weights Wp^ and the sum over representation 
labels p. This can be achieved by introducing another type of tensor = w, see figure 6. If we 
work with a cubic lattice these tensors are four-valent carrying as indices representation labels: 
(^■^)piP2P3P4 — Wp^5{pi, p2, P'i, Pi)- The second factor is equal to one if all representation labels 
in the argument coincide and zero otherwise. This tensor is connected by auxiliary edges to the 
four adjacent edge tensors ensuring that the sum over the representation labels involves the 
same representation label for every face. 

2.3 Relation between spin foam and spin net models 

Here we want to comment on an interesting relationship between spin net and spin foam models. 
Namely, one can understand spin nets as expectation values of (gauge symmetry breaking) 
observables inserted into the spin foam partition functions. We will illustrate this only for the 
simplest case: the spin foam represents BF theory and the spin net is of the form (2.1), i.e. the 
vertex projector is given by the Haar intertwiner. 
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Indeed, the partition function (2.1) can be rewritten in the form 



1 



z = 



{9e} e e 



1 



^En^^(5e)n^G(M' 



(2.32) 



oc 



\G\i' 



{9e} e / 



where / runs over a set of faces whose boundary edges and vertices generate the initial graph 
over which the spin net model is defined. Moreover, the two-complex made up by this set of 
faces must be simply connected for equation (2.32) to be valid (otherwise we have to amend 
the condition that holonomies along non-contractible loops should be trivial). The constant 
of propotionality comes from the normalization of the delta functions and can be absorbed by 
a redefinition of the weights We{ge) in the second line. The second line of equation (2.32) is 
the result of introducing the product of edge weights We{ge) in the partition function for BF 
theory {wj = 6g) on a simply connected two-complex. Similar observables have been considered 
in the context of 3D quantum gravity, namely for the Ponzano-Regge model [81, 82] and have 
been interpreted as Feynman diagram evaluations. That is, the edge- weights We{ge) can be 
understood as propagators and the spin net model as a Feynman diagram evaluation. 

3 Coarse graining methods 

Having introduced the models of interest in the previous section, we now turn to the challenge of 
implementing a coarse graining procedure. In this section, we will outline our general approach 
and discuss the conceptual issues that arise. The application of two approximation schemes in 
our context will then be discussed in the following sections. 

Although gravitational spin foams are usually defined on a general triangulation or two- 
complex we will here consider coarse graining on a regular lattice. This allows us to actually 
make explicit computations and to use methods from lattice gauge theory and condensed matter 
systems. One might object that using a regular lattice will introduce a background structure, 
spoiling background independence. There are, however, indications [25, 27, 83] that a restoration 
of diffeomorphism symmetry will be connected with a notion of triangulation or discretization 
independence, and hence the choice of a particular underlying lattice may not matter. Indeed, 
the universality phenomenon of statistical systems also suggests that the details of the chosen 
lattice might not matter for the questions we are interested in, i.e. a characterization of the 
possible phases of the models, or whether spin foams can avoid degenerate phases. Nevertheless 
one should study whether the results depend on the choice of lattice. 

The development of a scheme where order parameters or coupling strengths might be locally 
varying is another conceptual challenge [49, 50], which we will not address here. This would be 
appropriate for a random lattice or for situations with a very inhomogeneous dynamics. Again, 
we think that developing feasible coarse graining methods for spin foams and nets on a regular 
lattice is an indispensable first step. 

Furthermore, for gravitational systems, a regular underlying lattice can nevertheless repre- 
sent very inhomogeneous or irregular geometries. This is due to the fact that the dynamical 
variables are the geometrical data which also determine the lattice geometry and moreover the 
physical scale. In this sense a very fine lattice can carry very different (coarser) geometries 
[20, 84].^ This opens up the possibility that refining in lattice quantum gravity is even equiva- 
lent to summing over (a class of coarser) lattices [84] , thus eventually leading to a discretization 

®See also the universality result [85], which implies the simulation of irregular sublattices based on a regular 
underlying lattice. 
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independent theory at the fixed point [25, 27, 83]. Indeed, systems with diffeomorphism (like) 
symmetry might add interesting new insights to the theory of phase transitions. 

Another issue, which has to be addressed for any coarse graining or renormahzation scheme, 
is which space of models one considers the renormahzation flow to take place in [51]. Indeed, 
this question is not quite obvious for spin foam or spin net models as, for instance, spin net 
models mix aspects of edge models, where couplings are along edges, with vertex models, where 
couplings are based at the vertices, type of degrees of freedoms couple in a certain way might 
not be available for the non-trivial models. 

Furthermore, spin foam models are constructed to be as independent as possible from the 
underlying discretization. This translates into certain invariance properties of the amplitudes 
under a certain class of subdivisions, those which effect edges or plaquettes, but do not lead to 
an increase in the number of (non-trivial) vertices [70, 71, 86, 87]. Spin foam models, which 
are constructed using trivial face weights but non-trivial projectors, will satisfy this invariance. 
Invariance under this subclass of subdivisions can be seen as a first step towards a discretization 
independent model. One might ask whether it is possible to come up with a renormahzation 
scheme in which this form of the spin foam amplitudes and the invariance properties under face 
and edge subdivisions is preserved, see also [67]. However, this seems not to be very likely, as long 
as the amplitudes are not invariant under all kinds of subdivisions and hence renormahzation 
flow is trivial or at a fixed point. An intuitive reason is that the faces and edges of a coarse 
grained spin foam are 'effective' building blocks, containing a huge number of bare vertices, 
edges or plaquettes. That is, a change of the effective triangulation, even if it only involves a 
subdivision of edges and faces, will correspond to a more complicated change of the underlying 
'bare' triangulation. Indeed, we will find that under the renormahzation schemes presented here, 
the invariance property under edge and face subdivisions is not preserved. 

A general problem with real space renormahzation schemes of (higher than one-dimensional) 
statistical models is that an increasing number of non-local couplings appear with each blocking 
step.'' To keep the renormahzation flow in a space with a finite number of parameters, some 
sort of truncation or approximations scheme has to be used. We will consider two such schemes. 
The first one, the Migdal-Kadanoff [59, 60] scheme, is based on a truncation to local couplings. 
The derivation of this scheme is based on arguments that rely on the standard form of edge 
and plaquette models. Here an important question for future work would be to generalize this 
scheme to proper spin foam or spin net models, which constitute a rather generalized form of 
these models. 

Whereas the Migdal-Kadanoff scheme is based on a blocking of the group variables, the 
second scheme, based on tensor network renormahzation, is blocking the degrees of freedom 
encoded in the representation labels. This might be an advantage as the representation labels 
are related to metric degrees of freedom in the geometric spin foam models. Hence, this kind 
of blocking is much closer to the blocking used in the gravitation models in [25, 26], which is 
derived by geometrical arguments. Moreover we will see that this scheme allows direct access to 
the behavior of the vertex and edge projectors and makes use of the properties of representation 
theory encoded in the models. 

The tensor network renormahzation can, however, also be applied in the holonomy repre- 
sentation of the spin foam and spin net models, see also [90]. Which scheme to use might 
depend on which region one aims to explore, the high temperature (where the spin net/ spin 
foam representation is more appropriate) or the low temperature region (where the holonomy 

^Indeed, such non-local couplings are essential to regain diffeomorphism [26] or Lorentz symmetry [88, 89]. 
This at least applies to models with genuine field degrees of freedom, i.e. 4D gravity, but not for 3D gravity, which 
is a so-called topological field theory. Nevertheless, 3D gravity (which can be described by BF theory, that is 
the low temperature fixed point of lattice gauge theories) is an important test case for which the restoration of 
symmetries can be studied in a simpler setting. 
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representation might be more useful). Indeed, tensor networks are an extremely general tool, 
in which it is also possible to consider the emergence of different kinds of effective degrees of 
freedom, long range order as well as topological phases [91]. Another advantage is that tensor 
networks can handle non-positive weights [90] . This is an important point for gravitational spin 
foams where oscillating amplitudes appear, preventing the use of Monte Carlo simulations. 

The issue of non-local couplings in this renormalization scheme appears in the form of 
having more and more degrees of freedom. One has to choose a cutoff for this number. The 
advantage as compared to the Midgdal Kadanoff scheme is that the accuracy of this scheme can 
be systematically improved by increasing this cutoff. Furthermore, it is possible to study the 
dependence of the results on the choice of cutoff. 

On the other hand the tensor network renormalization scheme requires much more effort than 
the Migdal-Kadanoff scheme and offers less analytical control. This is of course related to the 
much bigger parameter space one is considering in the case of tensor networks. A crucial question 
for future research is how feasible the tensor network scheme will be for higher dimensional 
systems, as most work performed so far is for one and two-dimensional systems. Our application 
of the tensor network renormalization method is also restricted to two-dimensional spin net 
models. 

4 The Migdal-Kadanoff approximation 

Migdal-Kadanoff (MK) approximations are a simple tool to overcome a key difficulty in real- 
space renormalization: the introduction of non-local couplings in the renormalized action (which 
in the statistical physics language is the Hamiltonian) with each coarse graining step. 

When coarse graining the one-dimensional Ising model with nearest-neighbour interactions 
by decimation, the state sum factorizes. Integrating out the chosen spins gives a renormalized 
Hamiltonian of the same form as the one we had started with [60]. However, this no longer 
happens in higher dimensions: generically, the renormalized Hamiltonians feature new, longer- 
ranged interaction terms that render it impossible to continue with the procedure [59]. The 
central idea behind the MK approximation is to substitute this renormalized Hamiltonian by an 
approximate Hamiltonian which features the same interaction terms as the original Hamiltonian, 
thus making the renormalization transformation form-invariant. 

The original proposal [59] of Migdal amounts to neglecting certain terms in the state sum, 
while strengthening others. Applied to the two-dimensional Ising model, this translates to 
moving bonds (couplings along the edges) away from those nodes which are to be eliminated by 
decimation. In this way the valency of those particular nodes is reduced to two and the situation 
is thereby effectively reduced to the one-dimensional case (see figure 7). 

Kadanoff subsequently identified bond-moving as a special case of a general approximation 
scheme for which he derived a bound on the free energy [60]. Further justification for this 
approximation is derived from Monte Carlo simulations [52, 53]. In practice, MK approximations 
do fairly well at finding fixed points and phase transitions and not so well at determining the 
order of these transitions [54]. However, specific fixed points of Kosterlitz-Thouless type might 
not be reproduced by this approximation [92] . MK approximations are computationally efficient 
and comparably easy to implement but not refineable in a systematic way. 

4.1 Migdal-Kadanoff relations explored 

Let as consider a model of spin net type 
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(a) isotropic decimation 



(b) anisotropic decimation 



Figure 7: Derivations of MK relations for 2d spin net models involve moving bonds and sub- 
sequent integrations. The isotropic procedure leads to exponents 7 = 1,A = 2, whereas in 
the anisotropic case this is only true for the couplings in the vertical direction - the horizontal 
couplings feature 7 = 2, A = 1. 



based on the Abelian group TLq. Bond-moving then involves dropping some of the terms in (4.1) 
and replacing some of the others with 

w{g^,) w'^iQe') = ^ "^wik - j)w{j) XkiOe') (4.2) 
k \ j J 

leading to convolution of Fourier coefficients, the main feature of MK approximations. The 
subset of group variables corresponding to dropped terms can then be integrated out without 
generating non-local couplings. 

One can proceed similarly for lattice gauge theories or plaquette models of the form (2.19). 
For both kind of models different versions of bond-moving and different decimation schemes can 
be defined. In particular, in anisotropic methods the renormalized couplings will be different 
in different lattice directions, while in simpler isotropic methods the couplings will remain ho- 
mogeneous (see figure 7). We will only consider the isotropic methods here, which can also be 
made exact on so-called hierarchical lattices [93]. 

For the normalized weights Q{k) := w{k)/'w{0) this results in general recursion relations of 
the form [92] 

I E;:;«(-i)«(i) J 

Here, 7 and A are constants specific to the model and derivation in question. However, models 
with fixed 7 • A are qualitatively equivalent and share the same fixed point structure in the 

following sense: Two instances of (4.3), with exponents 7, A; 7, A resp. and 7 • A = 7 • A are 
related by 

Qnik) = [Qn{k)]^^'^ (4.4) 

given initial configurations that satisfy 

Qoik) = [Qo{k)r^^. (4.5) 

Typical values for the exponents are 7 = 1 and A = 2 for 'isotropic' 2D spin net and 4D lattice 
gauge models and 7 = 1 and A = 4 'isotropic' 3d lattice gauge theories. The corresponding fixed 
point equations are given by 

^J=0 / \i=o 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 O.C 0.8 1.0 



Q(l) = Q(3) Q(l) = Q(3) 

(a) A = 2, spin net 2D (b) A = 4, spin foam 3D 

Figure 8: Parameter space of Z4 symmetric models and their renormalization flow. The three 
regions correspond to configurations flowing to HTF (0, 0), LTF / BF (1, 1) and the 'cyclic fixed 
point' (0, 1). Other unstable fixed points (dots) and invariant submanifolds (thick lines) are also 
depicted. Note the differing behaviour of the cutoff-model (square) at (1,0) in both cases. 



The high and low temperature fixed points ( Q(0) = 1, Q{k > 0) = and Q{k) = 1 Vfc, resp.) 
are common solutions of eq. (4.6), independently of the exponents. Furthermore, for some fixed 
factor d of q, there is a class of fixed points given by 

I 1 k mod d = , , 

which derive this property from being invariant under d-fold cyclic permutations. 

More generally, for a fixed factor d of q, (symmetric and normalized) configurations with 
Q{k) = for k mod d ^ and those given by 

^ , I 1 k mod d = , , , 

Q{k) = { (4.8) 
I a else 

form invariant submanifolds of the parameter space. There are also more nontrivial fixed points, 
including the higher dimensional analogue of the Ising fixed point on the one-dimensional invari- 
ant line connecting the low and high temperature fixed points. For Z2, this point is predicted 
to be a nontrivial solution of 

Q(l + Q27)A _ (2Q7)A ^ 0. (4.9) 

In the 2D standard Ising model case (with exponents 7 = 1, A = 2) the solution is given by Q = 
wi/wq = 0.296 which corresponds via kT = Artanh((5)~^ to a temperature of kT = 3.282 (the 
exact solution is given by kTc = 2/log(H-\/2) ~ 2.269 [94]). Also note that MK approximations 
only predict this Ising- type fixed point for exponents with 7 • A > 1. Figure 8 illustrates the 
afore- mentioned features for 2D spin net and 3D spin foam models with group Z4. 

4.2 Restoration of BF symmetry in Abelian cutoff models 

Cutoff models are certain initial configurations for the MK renormalization group flow which 
are of particular interest to quantum gravity because they are derived from BF theory. The 
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question of interest here is whether the topological nature of BF theory, which is destroyed 
in the construction of the cutoff models, will be restored under the renormalization group flow. 
With a possible restoration of the BF phase, also the translation symmetries (2.14)-(2.31) would 
be restored, which in the 3D gravity models correspond to diffeomorphism symmetry. Whereas 
the BF phase — or low temperature fixed point (LTF) — in a 3D gravity setting represents flat 
space time, the high temperature fixed point (HTF), at which (5(0) is equal to unity and vanishes 
for all other labels k, corresponds to a geometrically degenerate phase, in which all geometric 
(length) observables have vanishing expectation value. However, whether the models flow back 
to BF or not depends very much on the initial conflguration in question. Instead of focussing on 
one model with one specific group, we here address the differences between different models for 
finite Abelian groups Zq with varying q, leaving other (non-Abelian) groups for further research. 

Let us first consider the case of 3D lattice gauge theory / spin foam models over Zq for varying 
q where the exponents in (4.3) are given by 7 = 1,A = 4. Here, different initial configurations 
parametrized by q and K converge quite fast (after 5 to 10 iterations) either to the HTF or to 
LTF, see figure 9b. In general, there are two competing effects encoded in the recursion relations 
(4.3): the convolution leads to a broadening of the function Q{k), whereas the exponents 7 and 
A lead to a dampening effect (as the Q{k) < 1). 

For the 3D gauge models most configurations fiow to the HTF, that is the dampening factor 
is quite strong. Regarding the question of restoration of the translation symmetries we see 
that this happens only in cases where the initial configuration is already quite close to the BF 
conflguration, which coincides with the LTF. 

These observations are in accordance with similar, but analytical work done in the generalized 
case of U{1) [92]. There it is found that for 3D U{1) gauge models all configurations satisfying 
certain conditions'^ flow to the HTF. 

These results have been extended to 3D lattice gauge models with non-Abelian compact 
groups U{N) and SU{N) by [95]. Hence this applies also to the 3D gravitational (Ponzano- 
Regge) model, which is based on SU{2) (assuming a change from a lattice based on tetrahedra 
to a cubic lattice does not matter) and we have to conclude that within the Migdal-Kadanoff 
approximation translation symmetry / diffeomorphism symmetry is not restored. Note, however, 
that for finite groups (i.e. finite q), there are some (even) cutoff configurations that do flow back 
to the LTF or BF phase. Here we might draw the conclusion that it is easier to restore a 
compact symmetry as compared to a non-compact one, as the translation symmetry is based 
on a compact parameter space for Zg and on a non-compact one for proper Lie groups. Here 
it would be interesting to see whether for some analogous modiflcations of the Tuarev-Viro 
models [33], which describe 3D quantum gravity with a cosmological constant and are based 
on quantum groups, such a restoration of the (here compact) translation symmetries occurs. 
See also a discussion of related issues in loop quantum gravity quantization of 3D gravity with 
a cosmological constant [34, 35] and a classical coarse graining treatment of the same system, 
where diffeomorphsim symmetry can be restored [25]. To study this question for the Tuarev- 
Viro models one would have to adjust the Migdal-Kadanoff method to quantum groups or to 
apply alternative renormalization schemes, such as the tensor network scheme described in the 
next section. 

For 4D 'isotropic' lattice gauge models, or equivalently for 2D spin net or edge models, 
the exponents in (4.3) are given by 7 = 1,A = 2. Hence the dampening effect is much less 
pronounced than for the 3D gauge models. Indeed, now most configurations flow back to the 
LTF or BF phase, see flgure 9a. An exception are the K = l,q > 8 configurations which flow 

^ These include a positivity requirement on the weight functions both in representation space and group 
space which is not satisfied for our cutoff models. However it would be satisfied for instance for a heat kernel 
regularization of for instance the Ponzano-Regge model. This would show that with the Migdal-Kadanoff method 
these regulated models would not fiow back to the Ponzano-Regge model, but to the high temperature fixed point. 
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Figure 9: Flow behaviour of different cutoff models, labelled by K and q. Markers half-filled 
at the top flow to BF/LTF (diagonal models, K = q/2), markers half- filled at the bottom flow 
to HTF (horizontal modesl with K = 0). In (a), unstable fixed points appear (pentagonal and 
hexagonal markers). 



to HTF, however only after a considerably large number of iterations (around 60). Nevertheless 
this can be interpreted as a phase boundary in the q — K diagram. 

For a given cutoff K > 1 the number of iterations necessary to converge to the LTF grows 
with q. (For K = 2 these are 8,15,60 iterations for q = 8,10,12 respectively). Moreover, for 
sufficiently large g, the simulations go through a long phase of only very small changes of the 
order 10~^ — 10~^, so that the iterations are almost constant. For K = 2 this appears starting 
with q > 14. For instance from around iteration 10 to iteration 100 the following configurations 
appear for the q = 14 simulation 

Q = (1, 0.81, 0.43, 0.15, 0.03, 0.005, 0.0005/6, 0.0001, 0.0005/6, 0.005, 0.03, 0.15, 0.43, 0.81) (4-10) 

and would be stable at least up to the number of digits displayed in (4.10). This configuration 
converges after 970 iterations to the low temperature fixed point.® Note however that this 
number of iterations corresponds to an extremely large lattice (in lattice units). 

This type of behavior is typical for fixed points with unstable directions but could also occur 
for quasi fixed points. To differentiate between these two cases, we considered a one-parameter 
deformation of the q = 14, K = 2 model, for which we changed the Q{1) = Q{2) = Q{12) = 
Q(13) values from unity to an arbitrary parameter < x < 1. Indeed, there is a phase boundary 
for X ~ 0.365 which leads to a non-trivial fixed point (we give only the first two non-vanishing 
digits) 

Q = (1, 0.72, 0.28, 0.057, 0.0062, 0.00035, 0.000011, 0.00000070, 0.000011, 0.00035, 0.0062, 0.057, 0.28, 0.72). (4.11) 

This fixed point is considerably different from the configuration (4.10) and also the initial con- 
figurations parameterized through x = 1 and x = 0.365 are quite different. Hence we see that 
a rather large portion of parameter space is dominated by the unstable fixed point. To obtain 
convergence to either LTF and HTF extremely large iteration numbers are necessary. In other 

®For higher q this convergence requires much more iterations, more than 36000 for q = 16. But the values 
(4.10) appearing through the stable phase would be very similar in the q = 14, g = 16 simulations. Indeed as will 
be explained in section 5 these two configurations can be considered to encode the same physical model. 
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words, the phase transition between LTF and HTF is very weak and the phase boundary not 
very pronounced. 

Indeed, in two-dimensional systems and in the Hmit of a continuous symmetry group, such as 
C/(l), one cannot expect the usual type of second order phase transition between the symmetry 
breaking phase (which here would be LTF) and the disordered phase (HTF). This is explained 
by the Mermin- Wagner (Coleman) theorem stating that continous symmetries cannot be spon- 
taneously broken at finite temperature [96-98] . Nevertheless, there are two phases, both in the 
two-dimensional system with global U{1) symmetry [99, 100] connected by a Kosterlitz-Thouless 
transition [101] and in the four-dimensional gauge system [102, 103]. This is a phase transition 
of infinite order. However, for the (isotropic) MK relations it was proven by Ito (again under 
certain assumptions on the initial configurations, including positivity of weights both in group 
and Fourier space), that this phase transition is not detected [92]. But also here extremely 
slow convergence (of all configurations towards HTF) has to be expected due to the existence of 
quasi fixed points. Our findings are explained by these considerations (although we use initial 
configurations not satisfying Ito's assumptions), that is by going to larger q we have to expect 
weaker and weaker phase transitions. Furthermore, the Ising type fixed point on the invariant 
line between LTF and HTF is drawn continuously towards HTF, indicating the absence of the 
transition in the limit of large q. For 4D lattice gauge theories with non-Abelian Lie groups 
it is speculated [92] that the MK relations provide a better approximations than for Abelian 
ones. The confinement conjecture states that in contrast to Abelian groups there is no phase 
transition for systems with non-Abelian Lie groups [103]. 

5 The tensor network renormalization scheme 

Here we will discuss a second coarse graining method and apply it to spin net models in 2D in 
the spin net representation. 

We have shown in section 2.1 that the spin net models can be easily brought into tensor 
network form. For these a number of real space renormalization techniques have been developed 
in recent years [61, 62, 104]. Moreover, tensor networks have became popular not only as a tool to 
formulate partition functions for 'classical statitistical models' but also to provide a variational 
ansatz for trial wave functions [73, 76, 105] for quantum statistical models. For a variational 
ansatz one has to find the expectation value of the Hamitonian with respect to the trial wave 
functions. The computational techniques [105] are similar to the tensor network renormalization 
group techniques. 

This is another point that motivates us to consider renormalization of spin net models, 
as structures very similar to spin nets, the so-called spin networks, appear in the canonical or 
Hilbert space formulation of spin foam models. Hence the renormalization techniques considered 
here could also be useful to find e.g. the physical wave function (ground state) via a variational 
ansatz. Indeed, the tensor network ansatz has also been developed to describe topological phases 
[62, 105], which often are represented by the so-called physical wave functions of i?F-theories, 
which are the starting point of spin foam quantization. Moreover, as we have also seen in 
section 2.1, tensor networks and spin networks are naturally related [74, 78]. More generally, 
tensor networks are a general tool for graphical calculus [78, 79], which becomes especially 
powerful for representation theoretic models such as spin foams and spin nets [77, 80]. 

The starting point of the tensor network renormalization (TNR) method is to write the 
partition function of a given model as a contraction of tensors associated to vertices of a graph 
(or lattice) 




(5.1) 



a,b,c,d... 
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The contraction of indices is along the edges of the graph, that is every edge of the graph carries 
one index. As these are summed over we can interpret the indices as the variables or degrees of 
freedom of the model. The dynamics is encoded in the choice of tensor. 

The tensor network itself, that is its underlying graph, can also be interpreted as a lattice in 
space or space-time. Choosing appropriate subsets of tensors and contracting all tensors inside 
each subset according to the connectivity given by the network will result in a network with 
fewer vertices and edges. Hence, this procedure corresponds to the blocking procedure in real 
space renormalization. Each subset results in a new 'effective' tensor T' describing an effective 
model. Notice, however, that in the two and higher dimensional case the effective tensors T' 
generically^" carry more indices than the original tensors T. For a regular lattice the number of 
indices associated to the effective tensors T' grows with the number of iteration steps. For the 
underlying graph it means that the valency of the effective vertices will grow - mostly in the 
form of having multiple edges between pairs of vertices. These multiple edges can be summarized 
into effective edges, which then carry indices with an exponentially growing range. 

Here is where one has to choose an approximation such that the indices run over a pre- 
chosen maximum number of values Dc- By increasing the approximation can be improved 
systematically. Ideally, this approximation should pick out only the relevant physics and neglect 
the irrelevant short distance fluctuations. The details of how this selection is implemented 
depend on the scheme, here we will follow a refined version of [61, 62]. 

The TNR method is very general as many different models can be written in tensor network 
form. That is in principle one can also flow between different models based on different kind 
of variables, models we obtain the same initial tensor network models for which also follow 
the same renormalization flow. On the other hand one loses a direct physical interpretation of 
the 'blocking procedure', i.e. how the effective/ blocked degrees of freedom are built up from 
the microscopic ones. This information can be supplemented by studying expectation values 
of (coarse grained) observables. Below we will introduce a method which allows to keep some 
physical interpretation of the indices associated to the effective tensors. This is however specific 
to the spin net/ spin foam representation, where the indices are group representation labels and 
for the gravitational spin foam models carry geometric information, such as lengths and area 
values. Hence we can argue that this method would correspond to blocking over these area 
or lengths variables. Such a blocking procedure was also employed in [25, 26] for the classical 
Regge model and has the advantage of a direct geometrical interpretation of the coarse graining 
procedure. 

Here we will consider the TNR method for a regular lattice but the principle is also applicable 
to generic graphs arising for instance from random triangulations or Feynamn graphs. In this 
case one needs, however, to find some suitable approximation scheme to prevent an exponentially 
growing index range of the effective tensors. 

5.1 Gaufi constraint preserving TNR method 

In this section we will shortly describe the TNR method following [61, 62] applied to 2D Abelian 
spin net models. We will, however, introduce a technique to keep the Gaufi constraints explicitly 
valid throughout the renormalization process, see also [74, 75]. The reason for doing this is that 
the Gaui3 constraints have an immediate geometrical information: In a given spin net (with 
oriented edges) consider any region such that its boundary cuts only through edges. Then only 
those configurations will contribute to the partition sum for which the sum of all ingoing indices 
is equal (modulo q) to the sum of all outgoing indices. This means that the Gaufi constraints 
should also hold at the effective vertices, which arise from blocking all the vertices in certain 

^"An exception are hierarchical lattices [93]. 
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(a) square lattice 



(b) splitting of vertices 



Figure 10: The square lattice with even (blue circles) and odd (green squares) vertices. 



regions. We will first review the method for a general 2D tensor network model based on a square 
lattice and afterwards specify to the case of spin net models and deal with the Gauf^ constraints. 

Consider a 2D tensor network based on a square lattice, so that the tensors T"*'^'^ are of rank 
four, see figure 10a. An obvious way to proceed would be to contract always four tensors along 
a square and to define in this way a new effective tensor which would now carry four double 
indices. 

However, to find a suitable approximation, that is a method to keep the index range constant, 
one proceeds differently. The first step is to decompose the tensors T into a product of two other 
tensors S. This is performed in two different ways according to the partition of vertices into odd 
and even ones. A vertex is even, respectively odd, if the sum of its lattice coordinates is even, 
respectively odd. 

For even vertices we decompose (see figure 10b) 

j.abcd ^ gab,igcd,i _ ^5 2) 

i 

Such a decomposition is always possible using a singular value decomposition (SVD) for the 
X d? matrix M'f''^'^ = T"'"^'^. Here d gives the range of the indices a,b, . . .. This gives 

^ab,cd ^ uf'^^iV^Y''^'^ (5.3) 
i=0 

with positive singular values Aj and unitary matrices U and V. We can then define 5^^'* = 
VXiUf^' and Sf'' = ^/AI(V^)^'^'^. 

Similarly for the odd vertices we decompose (see figure 10b) 

r,.abcd ^ Y sf^'sf^' (5.4) 

i 

where now one uses a SVD for the matrix M^''^'^ = T"-^'^'^, 

In a second step we contract four of the tensors S along the indices of type a, 6, ... to obtain 
the new tensor T'^^^\ now with indices i, j, . . . (see figure 11a), and arranged along a square 
lattice rotated by 45° (see figure lib) 

rp/ijkl ^ ^ gab,i gac,j gdc,k gdb,l 
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(a) contraction 



T' 



(b) coarse grained lattice 



Figure 11: (a) Contraction of the four S tensor to the new T' tensor, {h) The coarse grained 
lattice. 



If we keep the range of i as in equation (5.3) the index range of the tensors T would grow 
exponentially with the number of iterations. This is where the key approximation step comes 
in, namely to consider only the Dc largest singular values in the decomposition (5.3). This 
approximation is justified as the partition function is a trace over the tensors, thus involving the 
sum over the singular values. The validity of the approximation can be checked by comparing 
the values of the neglected singular values against the largest singular values in the SVD [61]. 
One can choose a rescaling after each iteration step such that this largest singular value is equal 
to one. Implementing the cutoff Dc in the number of singular values in the decomposition (5.3) 
we will obtain a flow in the space of tensors of rank four with a constant index range given by 

The SVD does not only serve as an approximation method but leads also to a field redefi- 
nition. Here the field variables are given by the indices over which the tensors are contracted. 
In the SVD these tensors are linearly transformed, which also induces a transformation on the 
fields. The transformations aim at an efficient representation of the partition sum, i.e. involving 
a minimal range of indices or equivalently minimal number or range of variables. The SVD is 
not unique, in particular for degenerate singular values one can add rotation matrices acting on 
the eigenspaces associated to the degenerate singular values. For the spin net models we will, 
however, modify the method so that the kind of field redefinitions that can occur are restricted. 
This is related to preserving the GauB constraints throughout the renormalization process. It 
also has the advantage that the indices keep their original physical interpretation, which for the 
gravitational models are related to geometrical quantities. 

Let us now specify to the Abelian spin net models in the spin net representation. Here it 
is convenient to introduce an orientation for the edges: for the square lattice we will choose all 
horizontal edges to point to the right and all vertical edges to point upwards. In this case the 
initial tensor T is of the form (the indices are anti-clockwise cyclically ordered starting from the 
leg pointing to the right, as in figure 10b) 

T''^^'^ = u{a)u{h)u{c)u{d) (5(«)(a + h-c 

where n(-) = w{-) in the notation of (2.16) and a = 0, . . . , (7 
delta function factor signifies the GauB constraints. Because of this factor the matrices Mi and 
M2 can be brought into block diagonal form, namely for M'^^''^'^ we have the condition that 

a + b = c + d=:i mod q (5-7) 

for non-vanishing entries, whereas for M^'"''^ we have that 

b-c = d-a=:j (5.8) 



d) 



(5.6) 



1 for a model based on Zg. The 
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Here the indices i,j label the non- vanishing blocks. 

In the first iteration step the decomposition into the tensors S can be obtained exactly and 
involves at most q non-vanishing singular values 

Sf' = u{a)u{b) (a + b-i), Sf'' = u{c)u{d) J^^^ (c + d-i), 

S^'^ = u{c)u{b) 6^i'> (b-c-j), Sf'^ = u{a)u{d) 6^'^^ {d-a-j) . (5.9) 

We assume that Dc > q (or alternatively that Dc is bigger than the number of non-vanishing 
u{a)), so that no approximation is necessary at the first iteration step. Note that at least 
Dc = q is necessary to flow to the low temperature fixed point, where u{a) = 1 for a = 
0, . . . ,q — 1. (As one can check the corresponding tensor is a fixed point also for the tensor 
network renormalization flow.) 

The contraction of four tensors S along the four edges of the square would involve four sums. 
Due to the (four) delta functions in the sum this, however, reduces to one summation. There 
will be one delta function 6^''\i + j — k — I) left, which amounts to the GauB constraint for the 
effective tensor T'*-^'^': 

T''^''^ = 6^''\i + j-k-l)J2u\c)u\i-c)u\j + c)u^{k-j-c) . (5.10) 

c 

This new tensor will in general not be of the factorizing form (5.6) anymore, so generically 
the decomposition into tensors S will now involve non-vanishing singular values. That is in 
case Dc < (f' , the approximation sets in. Nevertheless, the block diagonal form of all tensors 
and matrices involved can be kept through the following iterations. 

At a general iteration step we will work with a tensor T"'"^'^ with double indices a = (ja, ma) 
where ja = 0, . . . q — 1. As will be explained below the range of ttt-q depends on the SVD in the 
previous iteration step. These tensors will satisfy the GauB constraints 

']-{ja,ma)(h,mMjc,m,)(ja,ma) ^ (^-^ _^ _ _ j^-^ (^^ ^^^ 

at all iteration steps. 

Similarly as before we can deflne for the even vertices the matrix ^p"'™'")(j'*"™'6).(ic,mc)(id,md) _ 
j'{ja,ma){jb,mt){jc,mc)Ud,ina) _ GauB Constraint (5.11) this matrix can be brought into 

block diagonal form, as the non-vanishing entries must verify ja + jb = jc + jd ='■ i- We will 
denote these blocks by Mi(i). Then the singular value decomposition can be applied on the 
single blocks Mi(i) for i = 0, . . . ,q — 1, so that 

(^]\^^(^i'-^'^{ja,ma)ijb,mt),ijc,mc)ijd,m,d) _ ^ ■) (ia,ma)(jb,m6),(j,mO (^y^ ('^■jt ^(jc,mc)(jd,m ^ (5.12) 

rrii 

This yields of course the same singular values as for the entire matrix. Apart from providing 
a faster algorithm [74, 75] , the numerical implementation leads also to more stable results for the 
following reason. Generically one encounters the case of having singular values with multiplicities 
higher than one. In this case the singular value decomposition is not unique as one can perform 
rotations among the basis vectors associated to a given singular value. Here, keeping the block 
structure explicit prevents a mixing between different blocks induced by these rotations. 

To implement the approximation we now have to select the Dc largest singular values among 
the singular values of each block. The number Ni{i) of singular values selected from the block 
i determines the range of the second index rrii in the double index {i,mi). That is this number 
Ni(i) can take values between zero (no singular value selected) and Dc (all selected singular 
values come from one and the same block). Note that for the flrst iteration step Ni(i) = 1 (in 
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case all the u{a) are non-vanishing). We define the S matrices by 



gija,ma){jb,mb),{i,mi) ^ yO^((7]^(^))0'<""*<')0'i'>"*i')'(«>'"<) 

g{3a,ma){jb,mi,),{i,mi) ^ (■^)t^0c,"ic)0d,"id), _ (5.13) 

Note that we have i = ja +jb and i = jc + jd (mod q) for the non-vanishing entries of and 5*2 
respectively. For the matrices (M2(i))*^-''='"*=^^-''""*''^'^-''""*''^^^'*'"*<^^ at the odd vertices we proceed 
similarly, which will result in matrices and For the non- 

vanishing entries of these matrices we have jh — jc = i and jd — ja = i (mod q) respectively. The 
range of the indices r?ij is now determined by the number N2{i) of singular values selected from 
the block i of the matrix M2. 

The contraction of the four S-matrices along the square now results in a sum 

j,/{i,mi)(j,mj){k,mk){l,mi) _ ^(c,mc)(i-c,mi_c),(«,m»)^(c,mc)(i+c,mj_,_c),(i,mj) 

C 17lc,THi — c,tnj-^c^^k—j — c 

g{k-j-c,mk:-j-c)(j+c,mj+c),{k,mk) g{k-j-c,mk-j-c){i-c,mi-c),{l,mi) 

~ 5^'i\i + j - k - I) . (5.14) 

Here the range of mc, nik-j-c is determined by Ni{c), Ni{k — j — c) from the previous iteration 
step whereas rrii-c, rrijj^c is determined by N2{i — c),N2{j + c) respectively, also from the previous 
iteration step. Accordingly, the range of mi,mk is determined by Ni{i), Ni{k) of the iteration 
step under consideration (that is by the last SVD) and mj,mi by N2{j), N2{1) respectively. 

This altered algorithm has not only the advantage of keeping the Gauf5 constraint explicit 
but is also keeping any physical interpretation that might be attached to the representation 
labels / indices ja- For the gravitational spin foam models these would carry information on 
lengths and area variables - and the procedure here would correspond to a blocking where the 
microscopic geometrical variables arc basically added to obtain the coarse grained variables. 
The GauB constraints arise because of reasons rooted in representation theory: namely that 
for Abelian groups the tensor product of representations k and k' leads to a representation 
k + k' . Indeed the tensors T and S can just be seen as intertwining maps between (tensor 
products of) representation spaces, and the first index ja of the double index (ja, ma) encodes 
the representation carried by the associated edge. 

The same arguments based on representation theory apply for the non-Abelian spin net 
models. Also there the tensors T and S are intertwining maps between representation spaces 
leading to matrices of block diagonal form. For the abstract spin net models the initial tensor 
T is basically determined by the choice of projector in (2.9), which restricts the intertwining 
map from one of maximal rank to one of some smaller rank. Here it will be interesting to study 
whether any of the projector properties are preserved under coarse graining. 

Tensor network methods thus have the potential to give direct insight in the behavior of the 
projectors under coarse graining, which are the key dynamical entities in spin foam models. 
Moreover, the kind of coarse graining in the spin net or spin foam representation uses the 
representation theory underlying these models and furthermore corresponds to a geometrically 
natural coarse graining. 



5.2 Equivalence of models 

A crucial point in every renormalization method is the question which space of models one is 
considering, that is in which space the renormalization flow is taking place. Often the coarse 
graining process leads to models outside this space and usually some approximation method is 
employed to project the models back into the chosen space. For instance, in the Migdal-Kadanoff 
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approach one considers models with local interactions, that is non-local terms have either to be 
neglected or replaced by appropriately chosen local interaction terms. In particular for a Zg 
gauge model with local (single plaquette) interactions one always stays in the form of the 
gauge model (with single plaquette interactions). 

In contrast the space of models in the TNR method is defined by the chosen cutoff Dc on the 
index range of the tensors T. Hence different models, for instance spin net models with different 
Zq groups, can be considered in the same space. Indeed, for certain initial conditions the initial 
tensors T'"'"^'^ might actually define the same tensor network models. 

Consider, for instance, the Abelian cutoff models for different (even) q but for the same 
cutoff parameter K. That is we deal with the initial tensor (5.6) where 



with k running from — (| — l)to| (we consider even q). 

As T"^"^*^ = u{a)u{b)u{c)u{d)5^'^\a+b— c— d) the corresponding (symmetric) matrices M^''^'^ 
and Mg^'""^ will have zero rows and columns if these include an index a with u{a) = 0. In the 
SVD these rows/columns can just be neglected as it leads to a vanishing singular value. After 
neglecting zero rows and columns in the matrices arising from different models, these matrices 
might however coincide (with appropriate matching of indices) and in this sense define the same 
TNW model. 

Indeed, one can check that in the case of the models (5.15) this occurs for a fixed cutoff K 
but for varying q as long as g > AK + 2. In this case the first singular value decomposition for 
the matrices Mi , M2 give for both matrices the same AK + 1 non-vanishing singular values 

^max) ^max ~ Ij ^max ~ 1) ^max ~ 2, \max ~ 2, . . . , 1, 1 , where Xmax — 2i^ + 1 . (5.16) 

As long as the cutoff in singular values is chosen such that Dc > 4K+1 the first decomposition 
(5.2) of the T matrices into S matrices is exact. 

In particular this means that Abelian spin net models with q > 4K + 2 and fixed K but 
different q should go through the same renormalization sequence and hence should also end in 
the same fixed point. We will see in section 5.4 that this holds almost always in the numerical 
simulations. There will however be also examples where the fixed points depend on q. In these 
cases the simulations approach an unstable fixed point and the difference in the simulations 
appear only for large iteration numbers (more than 100 iterations). Hence, the difference can 
be explained by numerical instabilities and the fact that our method depends on the parameter 
q in order to keep the block structure of the T matrices explicit: q defines the number of blocks 
and hence the kind of possible field redefinitions which underly the method. 

Disregarding this point the equivalence between models should also hold if we send q ^ oo, 
that is consider U{1). In general, we see that the TNR method might also be applied to Lie 
groups (which would lead to infinite dimensional matrices) as long as we consider initial data 
such that the initial matrices reduce to finite dimensional ones due to either the appearance of 
zero singular values, or singular values which are sufficiently small. This equivalence between 
models appears only approximately in the Migdal-Kadanoff method, as there the number of 
parameters on which the renormalization fiow acts is fixed by q, or more generally the size of 
the group, on which the model is based. 

5.3 Structure of fixed points 

In the Migdal-Kadanoff scheme we considered a renormalization flow within a space described 
by q (or {q — 1) after normalization) parameters for spin nets based on the group Zq. We are 




1, for \k\ < K 
0, for |A;| > K 



(5.15) 
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much more flexible with the TNR method, where the number of parameters is determined by 
the chosen cutoff Dc on the number of singular values. Hence there are potentially many more 
fixed points. Indeed the fixed point structure is considerably more complicated, in particular for 
Dc » q, as we will exemplify below with the Ising model, q = 2. Note that (with the exception 
of Dc = 2, which can be treated analytically) we will only discuss fixed points which we found as 
a result of the renormalization process, i.e. by flowing to these flxed points. That is we will miss 
most of the unstable flxed points, which are those with repellent directions and would require a 
flne tuning of parameters to flow into. 

The main feature of the TNR method - with the approximation based on the singular value 
decomposition - is the appearance of non-isolated flxed points. These are argued [62, 104] to be 
due to short scale degrees of freedom which are not averaged out by the approximation method 
employed here. In [62] different forms of additional approximation steps (termed entanglement 
flltering) are suggested, that apply once the flow reaches the non-isolated flxed points. We will 
not consider these additional steps here and just describe the type of flxed point encountered 
in the original TNR method. For future work one should, however, study this issue in more 
detail. In particular, one could beneflt from a detailed comparison between approximations of 
Migdal-Kadanoff type and approximations based on the TNR method and furthermore from an 
analysis of relevant and irrelevant directions for the linearized renormalization flow around this 
kind of flxed points. 

In the following we will describe the results of the TNR method for the Ising model in the 
spin net representation (associated to the high temperature expansion) for the choice of different 
cutoffs Dc, see also [62] for a treatment of the Ising model with the TNR method improved by 
entanglement flltering. 

The smallest cutoff which would also allow a representation of the low temperature flxed 
point (LTF) n(0) = n(l) = 1 is Dc = 2. We start with conflgurations u{0) = l,u{l) = x. For 
X > 0.60858 these flow to the LTF represented by u{0) = u{l) = 1. Accordingly, the singular 
values Xi for this flxed point are Aq = 1, A = 1 with every block j = 0, 1 contributing one singular 
value. (We normalize the tensors in every step such that the largest singular value is equal to 
one.) 

For < X < 0.60857 the conflgurations flow towards the high temperature flxed point 
(HTF) represented by u{0) = l,u(l) = 0. This corresponds to having only one non-vanishing 
singular value Aq = 1- The transition point at u{l) = 0.6085 corresponds to a phase transition 
temperature of kTc ~ 2.572. This is a better approximation to the exact result of kTc ~ 2.269 
than the isotropic Migdal-Kadanoff result of kTc = 3.282 in section 4.1. 

In this case, Dc = 2, one can analytically compute the flow in some suitable parameters, for 
instance the tensor components T"'^'^'^. This allows to flnd a further (unstable) flxed point, given 
by the matrix '"^ = T"^'^'^ (the numbering of the rows and columns is ab = 00, 11, 01, 10) 



Ml-. 



l + s\ s?(l + si) 

■\{l + sl) 2s\ 

2sfs2 sl{l + sl) 

V s?(l + si) 2sls2 j 



,with si = 0.592902, si = 0.43105. (5.17) 



Note that the tensor describing this flxed point is no longer of the original form (5.6) as this 
would require S2 = sf. Nevertheless it leads to the phase transition between the low and high 
temperature regime. 

Let us now turn to the case Dc = 4. In this case conflgurations with x > 0.6466 will still 
flow to the LTF characterized by two non-vanishing singular values Aq = Ai = 1, one from every 
block. For conflgurations with x < 0.6465 we encounter however a new type of non-isolated 
flxed points, so-called Corner Double Line tensors [62, 104]. 
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Figure 12: Tensor with Corner Double Line structure 



'Double Line' indicates that we have to deal with pairs A = (a, a') of indices. The tensors 
are defined by associating a matrix C""^ to each corner of the four-valent vertex, see figure 12, 

rp{aa'){bb'){cc'){dd') ^a' b ^b' c' Qcd' Qda j^g'j 

Such tensors are fixed points of the renormalization flow (for Dc sufficiently large). The decom- 
position (5.2) of such a tensor is determined by the singular value decomposition of C"^. Namely 
if 

= ^u''^^ mv^'^ (5.19) 



we can write 



^(aa')m{cc')idd') ^ ^a'b u'' U^''' lUVl'V^' V^''' j C'^' . (5.20) 

Hence if the SVD of C"'' includes n non-vanishing singular values rjj, we will obtain n? non- 
vanishing singular values of the form ijjrn/ for the matrices obtained from T. 

The fixed point we encounter in the high temperature region for Dc = 4 is of CDL type with 
corner matrix 

c=(; °) . (5.21) 

That is in the SVD for the matrices associated to the tensor T we encounter four non-vanishing 
singular values Aj = l,y, y, y^. Here y ranges from zero and seems to get arbitrarily close to 
1 for n(l) = X reaching the transition point at 0.6496562.... Note that here we encounter a 
continuous fixed point family of 'high temperature type' ranging from having only one non- 
vanishing singular value equal to one to having four non- vanishing singular values, with three 
of these almost equal to one. If we take the transition between this fixed point family and the 
LTF (with two non-vanishing singular values equal to one) as the phase transition we obtain 
a critial temperature of kT^ = 2.221. This again is a better approximation to the exact phase 
transition temperature than the Dc. = 2 result. 

The appearance of these non-isolated fixed points is interpreted [62, 104] to be due to 
short range entanglement, which is not filtered out by the renormalization flow. Hence certain 
microscopic details of the models are remembered and the usual universality property of phase 
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transitions does not apply. One can nevertheless argue that the correspondence between phases 
and fixed points does hold: here models in a certain phase flow to a certain type of fixed points. 
Indeed we will see that for = 9 two different types of non-isolated fixed points appear, 
representing the low and high temperature regime respectively. In [62] additional approximation 
steps are suggested, designed to filter out this short range entanglement. We will leave the 
investigation of these and other additional approximation steps for future work. 

The appearance of the CDL type fixed points suggest that cutoffs = with n a natural 
number, might lead to a fast convergence, as these accommodate the v? non-vanishing singular 
values of a CDL type fixed points. We therefore will also discuss the case Dc = 9. Here, both 
the low temperature phase and the high temperature phase are described by non-isolated fixed 
points. For the initial value u{l) = x < 0.64293 we flow to fixed points, described by nine 
non-vanishing singular values of the form 

Aj = 1, y, y, y^, z, z, yz, yz, . (5.22) 

(For X = 0.64293 we have y = 0.7588 and z = 0.4200.) This suggest a fixed point of CDL 
form based on a 3 x 3 matrix C with singular values rji = l,y,z. This family of fixed points is 
continuously connected to the proper high temperature fixed point, where y = z = 0. 

For m(1) = X > 0.64295, that is in the low temperature phase, we also flow to a non-trivial 
flxed point, this time characterized by eight non-vanishing singular values of type 

Ai = 1,1,2/, ?/,?/,y,y^y^ ■ (5-23) 

(For X = 0.64295 we obtain y = 0.6434. For x ^ 1 the parameter y approaches and we obain 
the proper low temperature flxed point.) Hence this tensor cannot be of CDL type. Indeed it 
turns out that it is of the form 

T = TcDL Tltf (5.24) 

where Tcdl is a tensor based on the corner matrix (5.21) and Tltf is the tensor associated to 
the low temperature fixed point, that is of the form (5.6) with u{0) = n(l) = 1. Therefore eight 
non-vanishing singular values appear, as products of the four singular values A = l,y,y,y'^ for 
the CDL tensor and the two singular values A' = 1,1 for the LTF tensor. The appearance of 
this type of fixed point tensors of product form was conjectured in [62]. We will encounter more 
fixed points of this form in the next subsection. 

The transition between the high temperature family of fixed points and low temperature 
family of fixed points corresponds to a critical temperature of kTc w 2.274, which again approx- 
imates the exact result kT^ ~ 2.269 better than the = 4 result {kTc ~ 2.221). 

5.4 Analysis of Abelian cutoff models 

We will now discuss the renormalization behavior of Abelian cutoff models to compare it to the 
results obtained with the Migdal-Kadanoff method in section 4.2. We will consider the same 
kind of configurations as in section 4.2, parametrized by the size of the group q and by the 
cutoffii K. 

Now, as was pointed out in section 5.2 configurations with the same cutoff parameter K 
but different q might encode the same tensor network model. To utilize this we will choose 
the same cutoff Dc {Dc = 16, 25 and equal to 32 for some models) for different q, so that the 
renormalization flow should also be the same in these cases. This is different from the MK 
method, where the number of parameters always depends on q. 

^^Let us recall that K refers to the cutoff introduced in the space of group representations, that actually defines 
the model, while Dc is the maximum number of singular values considered when approximating the SVD. 
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Figure 13: Flow behaviour of different cutoff models, labelled by K and q, compare figure 
9. Markers half-filled at the top flow to LTF/LTFxCDL (diagonal models, K = q/2), empty 
hexagonal markers indicate a quasi fixed point or oscillating behaviour. Connected markers 
illustrate equivalent models in the sense of section 5.2. 



We have seen in section 4.2 that the MK renormalization fiow for the 2D spin net models 
was much more involved than for the 3D spin foam models: For sufficiently large groups 
configurations would undergo a behavior determined by unstable fixed points and rather weak 
phase transitions leading to very long convergence times. Therefore we have to expect similar 
properties to appear in the tensor network renormalization. Furthermore, as is usually the case 
if configurations approach a region around phase transitions, the approximation implemented 
by the cuttoff Dc might get less and less accurate. Indeed, the values of the neglected singular 
values might become comparably large (of order 10^^ of the largest singular value, which will 
always be normalized to one), even for the quite large cutoff Dc = 32. The TNR method, 
however, offers the possibility to increase the cutoff and thereby study the influence of the cutoff 
on the results obtained. 

Another issue that appears during a number of simulations are conflgurations for which the 
symmetry k — >• —k (equivalent to reversing all the edges) is broken. The reason for this is that 
singular values usually appear with a two-fold degeneracy, namely one from a fc-block, the other 
from a (— /i:)-block. The singular values from the block with A; = and k = q/2 are an exception 
to this rule. Also, quite often singular values appear with even higher degeneracy. Now, taking 
only a flxed number of singular values into account, one will quite often dismiss one singular 
value of a degeneracy pair and in this way break the edge reversing symmetry. As the number 
of A; = singular values, which are taken into account by the cutoff is not known a priori, and 
might even change during the iterations, it might be quite difficult to find a cutoff where this 
issue does not appear. 

In regions 'near' phase transitions, that is where convergence takes very long, these non- 
symmetric modes typically lead to unstable, oscillating behavior. This also signifies that the 
infiuence of the cutoff is not negligible. In the Migdal-Kadanoff method non-symmetric con- 
figurations may only appear due to numerical inaccuracies and it is quite straightforward to 
implement a symmetrization after each iteration step. A similar procedure for the TNR method 
(which is however not as straightforward to implement as for the MK method) would probably 
be very helpful to obtain more reliable results as well. Here we will interpret an oscillating 
behavior over large iteration times as indicating the presence of an unstable or quasi fixed point. 

We can broadly summarize the behavior we encountered in the simulations of the Abelian 



31 



cutoff models into the following classes, wich we also illustrate in figure 13: 

• The models flow quite fast (typically during 10 to 30 iterations) to the low temperature 
fixed point or to a fixed point of type Tltf ® Tcdl as described in section 5.3. This 
happens for cutoff models which are not independent of q. The low temperature fixed 
points come with q non-vanishing singular values equal to unity. The more complicated 
fixed points with CDL structure appear, in particular, for the higher cutoff = 25. 
The singular values associated to the CDL structures are however quite small 10~^ 
to 0.2), so that we can safely interpret these fixed points as of low temperature type. It 
also happens that initially the configurations fiow to a T^tf 'Si Tcdl fixed point with a 
higher number of non-vanishing singular values than provided by the chosen cutoff. In 
this case the configurations slowly keep changing such that the singular values associated 
to the CDL structure decrease. With the exception of two models {K = 3, g = 12 and 
K = 4:,q = 16) such a behavior appearing in the D(. = 16 simulation would be confirmed 
in the Dc = 25 treatment. 

• An interesting special case is K = 1 and q> 6 (which describes the same tensor network 
model independent of q as long as g > 6). For Dc = 16 the fixed point (reached after 
about 60 iterations) is of the form Tltf "X" Tcdl, but the Tltf factor comes with only 
four singular values equal to unity. The CDL structure leads to 16 non-vanishing singular 
values, where apart from the four singular values equal to one, the next eight ones are 
equal to 0.86 and the final four equal to 0.73. 

The Dc = 25 results are slightly different: due to the appearance of non-symmetric config- 
urations as described above, the configurations are oscillating for more than 80 iterations. 
Then the results do actually start to depend on q: for g = 6 a Tltf ® Tcdl fixed point 
is reached with six singular values equal to one (so it is really a g = 6 low temperature 
fixed point), 12 singular values equal to 0.81 and six equal to 0.655. Whereas for g = 8 a 
fixed point is approached (around 200 iterations) with eight singular values equal to one, 
another sixteen are around 0.3, plus eight singular values around 0.12. This gives more 
than 25 non-vanishing singular values, so as described above the configuration is slowly 
changing, decreasing the CDL singular values. Note however that due to the high iteration 
numbers involved and the appearance of 'non-symmetric' configurations the results should 
be taken with some care. 

• For the Dc =16 simulations there are two examples K = 2 and q = 10 or higher, and 
K = 2i and gr = 14 or higher, which show stable behavior and seem to approach a (non- 
trivial) fixed point for long iteration times (approx. 20 iterations for K = 2 and approx. 
100 iterations for K = ^i). The change in the singular values during this approach can 
get very small, with instances where the change in the singular values are of the order 
10"'^. There is a large number of non-vanishing singular values indicating a complex 
(quasi/unstable) fixed point. After this stable phase configurations become unstable and 
'non-symmetric' showing a slightly oscillating behavior (changes of the order 10~'^ — 10~^ 
from one to the next iteration). For K = ?> the configurations converge to g-dependent 
low temperature fixed points (at iterations 230 for g = 14 and 180 for q = 16), for iiT = 2 
the behavior remains unstable for very long iteration times. This behavior is confirmed 
for K = 2 with the Dc = 25 and Dc = 32 simulations (where now the stable phase 
is up to 100 iterations). For K = 'i, the Dc = 25 and Dc = 32 simulations lead to 
non-symmetric configurations which show oscillating behaviour for long iteration times. 
However, the complicated structure of the largest q singular values agrees with the one of 
the simulations with smaller cutoff. 
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Note that examples flowing to a high temperature fixed point do not appear. The best 
candidate would be the K = 1 configurations (describing the same tensor network models 
starting with q = 6 and larger q), which however flow for Dc = 16 to a fixed point with four 
singular values equal to unity, plus more non-vanishing singular values due to the CDL structure. 

In the following table we list the details of our findings and compare them with the results 
obtained with the Migdal-Kadanoff method. 



K 


1 


MK 


TNR Dc = 16 


TNR Dc = 25 


1 


4 


LT(oB) 


LT(6)XCDL(7) 


LT(6) X CDL(>ioo) 


1 


6 


LT(i8) 


LT4(60)XCDL(63) 


Osc^LT6(42) X CDL(48) 


1 


8 


HT(59) 


LT4{60)XCDL(63) 
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Li (04) 
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1 /I 
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3 


8 
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3 


10 
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3 


12 
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14 
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16 
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LT(i4o) 




Osc 


3 
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16 
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Table 1: Summary of the renormalization flow of Abclian cutoff models. To a precision of 10~^°, 
the numbers in brackets give the iteration times (halfcd for TNR to make comparable) it takes 
to reach a certain behaviour: LT (low temperature fixed point, bold numbers denote the number 
of ones that appear explicit), HT (high temperature fixed point), QF (quasi fixed point), Osc 
(oscillating behaviour), CDL (corner double line), NULL (finite part that dies off). For a given 
K and growing q, the fine line indicates the beginning of physically equivalent models in the 
sense of section 5.2. 

The overall picture of the Migdal Kadanoff results is confirmed by the tensor network sim- 
ulations: Most configurations flow to the low temperature fixed point or to a low temperature 
fixed point embellished with a CDL structure. With growing q, cutoff models with sufficiently 
small K show long stable phases where a fixed point seems to be approached but then enter an 
unstable phase with (slightly) oscillating behaviour. For some cases these converge finally to a 
LTF fixed point, however this result has to be taken with some care, due to the oscillating phase 
in which 'non-symmetric' configurations appear. 

There is one important difference between the MK and the TNR results. As we explained 
in section 5.2 examples with q > 4K + 2 should encode the same physical models both exactly 
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and in the approximation provided by choosing the cutoff Dc- This property is inherent in 
the TNR method (and only appears to be violated for large iteration numbers due to the 
oscillating behavior). But in the MK method the truncation is determined by q and hence the 
renormalization flow is still different for q not much larger than qx '■= + 2. For much larger 
q the renormalization sequences will turn out to be almost the same, however this behavior sets 
in later than in the TNR examples. Moreover from [92] we have to expect that going to g — t- oo 
the configurations should flow to the high temperature fixed point. 

This is the reason why the quasi fixed point behavior has to set in earlier for the TNR 
method, as here the results starting with qx should in principle also hold in the limit g — )• oo. 
From this point of view it is interesting that we have not seen an example in the TNR method 
(especially not K = 1, as opposed to the MK K = l^q = ^ example) which would flow to the 
high temperature fixed point. Hence the TNR method might be able to detect two different 
phases for the U{1) theory. 

As we have also seen definite conclusions are very much hindered by the appearance of non- 
symmetric configurations, i.e. where the k~ blocks would differ from the — fe-blocks. These also 
appeared in the MK simulations for the examples which go through very long almost stable 
phases. In this case non-symmetric configurations only appeared due to numerical errors and 
this problem can be easily cured in the MK method by symmetrizing the Q{k) parameters after 
each coarse graining step. 

In the TNR method non-symmetric configurations not only appear due to numerical errors 
but are also caused by the cutoff, which might neglect one of a pair of degenerate singular values. 
This happens quite generically if the cutoff is increased (apart from the Dc = 25 simulations we 
also tried Dc = 32 and inbetween values) and it is increasingly difficult to find a value for Dc 
where this does not appear, say, for the first twenty iterations. 

This non-symmetric behavior does not matter so much for configurations which would con- 
verge fast to some stable fixed point but do cause long sequences of oscillating behavior for 
configurations, which we suspect would otherwise rather approach slowly an unstable fixed 
point. For these examples the ordered sequence of singular values at a given iteration decreases 
rather slowly. Even for Dc = 32 the neglected singular values can be of order 10~^ and an 
unsymmetric cutoff will have considerable influence on the overall behaviour. This problem is 
especially pronounced for the 2D models with larger g, as we there encounter rather weak phase 
transitions. 

For future work one should address this issue^^. There are different possibilities, one is to use 
a symmetric parametrization in the coarse graining procedure, i.e. only work with the A:-blocks 
where k < ^. Another option is to use an adaptive cutoff Dc for each iteration step, such that 
it avoids to cut between (k, —k) pairs. 

A third option would be to change the approximation scheme slightly by choosing a block 
dependent cutoff Dc{k). This would actually simplify the algorithm considerably. For sufficiently 
high Dc and Dc{k) these schemes should yield equivalent results, this has however to be tested. 

Note that non-symmetric configurations can however also appear due to numerical instabili- 
ties, so that one might have to implement some symmetrization procedure (if one does not work 
with a parametrization that only allows symmetric configurations). 

6 Discussion and outlook 

Extracting large scale physics from spin foam models is one of the most pressing issues for 
the field. Here we advocated the development and use of coarse graining and renormalization 

^■^It does however only appear for groups where the dual representation p* is not equivalent to the original 
representation p. For the rotation group SU{2) we rather have p* = p, but for (7(1) we have pi = p-k- 
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techniques so that hopefuUy a feasible method for 4D models can be found. To this end, we 
introduced a wide range of simplified models, spin foams with finite groups and moreover spin 
nets, which can be seen as a dimensionally reduced version of spin foams. This not only enables 
us to test and develop coarse graining methods but also to obtain some physical insights and to 
put forward conjectures about the dynamical behavior of the full models. 

This strategy allows us to adopt coarse graining methods from lattice gauge theory and con- 
densed matter system, here the Migdal-Kadanoff scheme and the tensor network renormaliza- 
tion method. The two methods have different drawbacks and advantages. The Migdal-Kadanoff 
scheme facilitates quick results (even a large number of iterations can take only seconds on a 
PC), so that an overview of the phase structure encoded in the model can easily be obtained. 
Here the main question is how this method can be generalized to non-Abelian spin foams with 
non-trivial projectors, which no longer fall into the class of standard lattice gauge theories. 

The tensor network renormalization method has the advantage of providing a systematic 
improvement on the accuracy of the results. The required effort is considerably larger (100 
Dc = 32 iterations may take several days on a PC). The method is however very general and 
in its version based on the spin net representation, allows direct access to the behavior of the 
(vertex) projector under coarse graining. Moreover, the blocking of variables is very natural if 
one takes into account the geometric meaning of the representation labels in the gravitational 
spin foam models. We presented an algorithm in which the GauB constraint are kept explicitly 
intact. 

The tensor network renormalization method is easily generalizable to models with non- 
Abelian groups, and indeed would nicely interact with the group theoretic content of these 
models, see also [74]. Here it will be very interesting to study how a non-trivial vertex or edge 
projector might change the phase structure as compared to the standard choice of the Haar 
projector. This will also facilitate a better understanding of the dynamics in the full gravitational 
models. To this end, a class of finite group models emulating the current gravitational EPRL 
models [8] is constructed in [106]. An important question for future research will be whether the 
degenerate phase, that is the high temperature fixed point can be avoided by selecting suitable 
projectors. This can already be studied for the 2D spin net models. Furthermore it has to be 
explored how the tensor network renormalization method can successfully be applied to three 
and four-dimensional spin net and spin foam models. 

An alternative to using the tensor network formalism to describe the partition function of 
a system would be to apply tensor network renormalization as a kind of improved mean field 
approach in a canonical quantization. This would change the (statistical) systems from D~ 
dimensional classical to {D — l)-dimensional quantum ones. In this case the tensor networks 
would provide an ansatz for the ground states of the models [73, 76, 105]. In gravity, instead 
of minimizing one Hamiltonian or energy functional, one rather has to deal with a number 
of constraints, which have to annihilate the so-called physical states. The master constraint 
[107-109] or the uniform discretization approach [110, 111] provide a framework where just one 
(master) constraint has to be minimized. 

In this work we have considered systems on a regular lattice, as this made explicit simulations 
feasible. Nevertheless one should consider generalizations of these methods to random lattices 
[49, 50]. Another question is how the phase structure found on a fixed lattice relates to phases 
in models, that include a sum over all possible lattices [112-115], in particular regarding the 
remarks in section 3. 

We see this work as a contribution towards a closer link between the quantum gravity and 
the statistical physics / condensed matter communities. For quantum gravity researchers, this 
offers the prospect of new concepts and numerical tools to study the large-scale physics of their 
models, whereas for people working in statistical physics it offers new models, new questions, a 
geometrical perspective and a rich set of mathematical tools behind it. 
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